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Abstract. A simple ordinary differential equation is derived governing 
the red-shifts of wave-fronts propagating through a non-stationary 
spherically symmetric space-time. Approach to an event horizon 
corresponds to approach to a fixed point; in general, the phase portrait 
of the equation illuminates the qualitative features of the geometry. In 
particular, the asymptotics of the red-shift as a horizon is approached, 
a critical ingredient of Hawking's prediction of radiation from black 
holes, are easily brought out. This asympotic behavior has elements in 
common with the universal behavior near phase transitions in statistical 
physics. The validity of the Unruh vacuum for the Hawking process 
can be understood in terms of this universality. The concept of surface 
gravity is extended to to non-stationary spherically symmetric black 
holes. Finally, it is shown that in the non-stationary case, Hawking's 
predicted flux of radiation from a black hole would be modified. 
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1. Introduction 

By definition, a black hole cannot be directly observed by someone in its exterior. 
For this reason, it is the geometry of the neighborhood of the horizon, as measured by 
later and later observers, which is of interest in the external physics of the holes. A 
basic element of this geometry is ray-tracing: looking in from far away towards where 
the hole will form (and neglecting opacity), what would one see? 

An estimate of just this, for rays sufficiently close to the horizon in the 
Schwarzschild (or more generally Kerr-Newman^) geometries, is a cornerstone of 
Hawking's (1974, 1975) prediction of thermal radiation from black holes. In the 
Schwarzschild case. Hawking argued that the asymptotic ray-tracing relation between 
surfaces of constant phase labeled by retarded time u in the future to those of constant 
phase labeled by advanced time v in the past would be 

t)(u) ~ Cexp — u/(4m) , (1) 

as tt — > +00, that is, as the horizon is approached. Here m is the mass in geometric 
units. It is v{u) which is the key function in the Hawking analysis; for example,the 
Hawking temperature comes out of the calculation as Th = —{l/2n)v/v. 

The present paper is concerned with the justification of the relation (1), and the 
derivation of similar relations for more general situations. It will be well to first explain 
what is known. 

Hawking's argument for the relation (1) was brief and perhaps elliptic; the taking 
into account of the effects of tracing the rays back through the collapsing matter was 
not really spelled out. The argument relied on specific properties of the Schwarzschild 
(or Kerr-Newman) solutions. A fuller treatment of the spherically symmetric case, 
assuming a static exterior, was given by Birrell and Davies (1982). Their argument 
was by computation in a particular coordinate system, and its relation to invariant 
geometry is not manifest. 

These arguments leave open two sorts of issues. First, one would like an invariant 
explanation of where the relation (1) comes from; and second, one would like to be 
able to treat a broader class of objects. For example, one would like to know what 
happens for a black hole accompanied by an accretion disk, or a spherically symmetric 
hole by a (non-static) nebula. What happens if the collapsing matter is not sharply 
bounded (so the event horizon may not be in a static exterior)?^ If we are concerned 
(as we are in the case of the Hawking process) with a star that has contracted to 
within atomic or even nuclear dimensions of its Schwarzschild radius, are we justified 

^ In what follows, it will be understood that the Reissner-Nordstrom and Kerr solutions are special 
cases of Kerr-Newman. 

^ There is a loose sense in which many astrophysical black holes are expected to be Kerr Newman 
("black hole uniqueness"). However, this is because for this class of holes one assumes that the hole 
approaches a stationary state and is isolated (that is, has only the electromagnetic stress-energy in its 
exterior). 
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in treating it as a classical bounded object? And what about bodies which are not 
quite black holes? What happens if a body collapses to r = 2m(l + e) and stays there? 
To investigate conceptual questions about the Hawking process, it is desirable to have 
an invariant way of understanding and treating these issues. 

The aims of the present paper are: (a) to bring out very clearly where the 
formula (1) comes from; (b) to show that it has a flavor of "universal" behavior 
in the sense of statistical physics; (c) to show that a similar result will hold without 
any assumption of asymptotic stationarity or vacuum exterior for the hole; and (d) to 
present a general framework for the study of similar questions. 

In the course of this analysis, we clarify some general features of the structure 
of black holes. We shall find a way of extending the concept of surface gravity 
to non-stationary spherically symmetric black holes. This new measure of surface 
gravity will be a global one, determined by red-shifts of light rays passing close to 
the horizon. While global, it is given by a well-defined integral formula in terms of 
the local geometry. This enables us to ask which details of the geometry the surface 
gravity is sensitive to: which contribute to the integral? 

We shall see that, as the hole is approached, all contributions from the regime 
prior to the formation of the event horizon are suppressed, and so the surface gravity 
is independent of the details of the formation of the hole, as one would expect. 

We shall also see that the surface gravity is essentially insensitive to the geometry 
of the horizon! This is in apparent confiict with the usual view that the surface gravity 
is determined by the gradient of the norm of a Killing vector at the horizon. There is no 
real conflict, however. The usual definition of the surface gravity is really a global one, 
since Killing's equation serves to evolve the Killing field throughout the space-time, 
and the Killing field must be normalized at infinity. Thus Killing's equation provides 
a sort of rigidity which connects data in one part of the space-time to data in another. 
This rigidity makes it impossible to use the Killing field to localize contributions to the 
surface gravity. Our formula, on the other hand, expresses it as an integral, and one 
can say which elements contribute more substantially than others. What the surface 
gravity turns out to depend on is how the neighborhood of the horizon is linked to 
future null infinity, not on what happens at the horizon itself. 

Of course, in the stationary case, the understanding of the surface gravity 
developed here, and the usual one, are simply alternative ways of regarding the same 
thing. In the non-stationary case, we have a new candidate definition for the surface 
gravity. One would not expect this to be constant. For example, if a black hole were 
enclosed in a shell of matter (or radius far larger than the gravitational radius of the 
hole), one would expect all physics within the shell to be red-shifted due to the shell's 
potential. Thus the surface gravity of the hole should appear smaller from outside the 
shell than if the shell were not there. If the shell were moved slowly ("adiabatically"), 
one would expect to a time-dependent surface gravity, with time-dependence given 
by the time-dependent red-shift due to the varying potential of the shell. 

Our formula for the surface gravity is well-defined in such circumstances, and 
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beyond the adiabatic approximation. It remains independent of the details of the 
formation of the hole, and (in the sense of the second paragraph previous) of the 
geometry of the horizon. 

The techniques we use are fairly elementary ones from differential geometry and 
ordinary differential equations. We shall show that the function v{u) is governed by a 
certain system of ordinary differential equations. (These equations are equivalent to 
some of the familiar optical equations.) Approach to an event horizon corresponds to 
approach to a certain kind of fixed point in the phase diagram; the asymptotics of the 
equations near the fixed point are what lead to the constancy of the surface gravity. 

Most of this paper is concerned with the spherically symmetric case. Section 2 
describes the universal character of the asymptotic behavior of the red-shift. Section 
3 derives the main equations governing the evolution of the red-shift, and Section 4 
discusses their consequences, in particular, the behavior of red-shifts near an event 
horizon. Section 5 contains comments on the non-spherically symmetric case. Section 
6 sketches the implications for Hawking's proposal in the case of non-stationary black 
holes. 

Conventions. The conventions are those of Penrose and Rindler (1984-6). The 

metric has signature H , and the curvature satisfies Va VfeV*^ — VbVaV'^ = Rab/V^- 

We use natural units. 



2. The Universality 

I will begin with a simple physical picture which brings out the "universal" 
character of the red-shift, without actually establishing its functional form. (Precisely 
what is meant by this will appear shortly.) This argument contains some restrictive 
assumptions, but its simplicity makes it worthwhile. 

Consider a spherically symmetric collapsing object in general relativity. We shall 
assume the object has a sharp boundary, or limb, and that exterior to this is vacuum. 
Then the exterior will be a portion of the Schwarzschild solution. Inside the collapsing 
object is another metric, the details of which will be unimportant. We will only need 
to assume that the interior metric joins to the exterior metric suitably smoothly. (A 
metric of class would be more than enough.) We shall assume a black hole forms, 
so that as t — > +oo the coordinate r{t) of the limb will satisfy r{t) — ^ 2m. 

Now trace backwards in time a radially-directed wave-packet. The wave packet in 
the distant future (near future null infinity, X"*"), has some characeristic wave-length 
with respect to the frame defined by the spherical symmetry at X+. We are going to 
follow it backwards in time, through the collapsing object, and out to the distant past 
(T~). We do this in two stages, which have different physical significances. 
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Figure 1. Penrose diagram 



for tlie splierically symmetric collapse of a massive body (shaded). The event horizon 
is Ji. A geodesic slightly before Ji is also shown, divided into two stages along which 
we consider the propagation of a wave packet. 

Stage I. This covers tracing the packet back from T"*" to the hmb of the coUapsing 
object, which is (say) at coordinate r(tpacket out)- If ''(^packet out) - 2m, then the 
wave packet is deep in the potential weU and has been very red-shifted. The 
precise measurement of the red-shift depends on the choice of the observer's frame at 
r(tpacket out)- However, relative to, say, observers freely falling into the hole, the red- 
shift is large, and becomes larger as later and later packets are used, since for these 
^(^packet out) 2m. As later and later packets are used, this portion of the trajectory 
corresponds to a divergent red-shift. 

Stage II. This covers tracing the packet back through the collapsing object, out 
the other side, and off to past null infinity X^. We would expect this propagation to. 
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to some extent, undo the red-shift of the journey from to the hmb. However, it 
cannot undo all of the red-shift, in the following sense. Throughout the second stage of 
the trip, there is no propagation which can give rise, even in the limit of later and later 
times, to an infinite red-shift. On this portion of the trip, the red-shift as the packet is 
propagated through the collapsing body remains finite, since this is only propagation 
over a finite non-singular region. Also the propagation from the collapsing body to 
X~ introduces no divergent red-shift, since on this segment of the trip, the radius of 
the body stays bounded away from 2m so no infinitely deep potential well is involved. 

The total red-shift is the product of the red-shifts for Stages I and II, and thus 
has the form of a part which diverges (as later and later packets are used) times a part 
which remains regular. The divergent part depends only on the propagation through 
the vacuum region exterior to the body, and so is independent of the details of the 
collapsing body, that is, universal for the class of spherical sharply-bounded bodies of 
given mass. 

We may extract the invariant information in this universality class as follows. 
Write the red-shift factor as v{u) = wiwu for the product of the two stages; then, 
since wu will tend to a constant, we see that the limiting behavior of (log?))' = v/v 
is independent of the trip in Stage II. It is also independent of the particular choice 
of Lorentz frame for measuring red-shifts near r = 2m, since changes of frame will 
only contribute smooth multiplicative factors with finite limits. At the present level 
of analysis, the limiting form of (logv)" could in principle be anything of the form 
f{u/m)/m, where u is the retarded time and / is a dimensionless function. Explicit 
computations show that — (log?))' — k = l/4m, the surface gravity of the black hole. 

While we have emphasized that there is an element of physics here in common with 
the universal behavior familiar in statistical mechanics, it would be wrong to think 
that the parallel extends, in any obvious way, very much further. The picture drawn 
here of the formation of a black hole differs from that of a conventional thermodynamic 
system fundamentally. First, we are concerned here with an essentially non-stationary 
process, since the time-dependence of the limb's trajectory is crucial. Second, the 
space-time is not a homogeneous extensive system. For example, one would not have 
a correlation length in the usual sense, although at a formal level one can treat the 
red-shift factor v as such. 

3. The Case of Spherical Symmetry 

In this section we derive the differential equations governing the evolution of the 
red-shift, in the case of spherically symmetric space-times. 

We shall assume (M, gab) is a spherically symmetric space-time, which is 
asymptotically fiat at past and future null infinity in a sense strong enough to 
guarantee the existence of X = U T~ with their usual differential-geometric 
structures. Choose Bondi systems (tt, 9, (p) on X+ and (v, 9, (p) on I~ respecting the 
spherical symmetry. Then u will be called the retarded time coordinate, and v the 
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advanced time coordinate. Fix too a null vector field n" tangent to the generators of 
and a null covector field orthogonal to the Bondi cuts at T+, with la'n"' = 1. The 
normalization n"Vatt = 1 is assumed. 

Caution. In most papers it is conventional to use for a tangent vector to X+ 
and for a tangent vector to X~. We shall not do this. The arguments here are global 
ones involving parallel transport between I~ and X"*", and it is most natural to take 
and defined at X+ and then extend them to X~ by transporting along certain 
paths. This gives the opposite-to-usual senses for the vectors on X~. 

Note. A rigorous treatment of the differential geometry at X would involve 
conformally rescaling the metric and also rescaling the components of tangent vectors 
so that their limits took finite values on X. This seemed overly formal for the relatively 
simple geometric issues to be considered here. Thus all calculations in this paper "at 
X' should be understood as done close to X, in the asymptotically flat regime. The 
sophisticated reader will have no difficulty recasting the results and arguments here 
in a more formal mold. 

3.1. Parallelism at Infinity 

One might hope that very far away from an isolated system one can use parallel 
transport to provide an essentially unambiguous way of identifying vectors at one point 
with vectors at another. In general, the situation is not quite so simple (because, while 
the gravitational field is locally weak, the regime far away from the source contains 
points very distant from each other). It turns out however that this is a reasonable 
hypothesis in the spherically symmetric case. (It is gravitational radiation which 
interferes with the use of aparallel transport to identify vectors at different points of 
X; in a spherically symmetric space-time, there is no gravitational radiation.) In other 
words, we assume that over X = X+ UX^ we have a natural way of identifying vectors 
at different points, which we call a parallelism. 

One important conseqence of this is that the red-shift of a radial null geodesic 
is well-defined. Suppose that we send in a wave packet with wave-vector If^ in the 
distant past (at X~), and it emerges as l^^^ at X+. In order to compare If^ and /"^^ 
we must have some way of identifying the tangent spaces at the point of emission on 
I~ and the point of arrival on X"*". It is this which the parallelism provides. Using 
the parallelism, we can regard /^^ and l^^^ as elements of the same space. Then the 
symmetry of the situation implies that that must be multiples of each other, and the 
red-shift is simply the constant of proportionality. 

It should be noted that the red-shift factor can be given directly. Take a-u — const 
surface near X"*" and trace it backwards until it emerges as v = v{u) = const near T~ . 
Then the function v = v{u) is the mapping of surfaces of constant phase, and v{u) is 
precisely the red-shift factor.^ 

^ While it may seem that v{u) gives a definition of the red-shift independent of the existence of a 
paralleUsm, this is not really so. For in order to interpret v{u) as a red-shift, one needs to know how to 
relate the clocks near X+ (relative to which the coordinate u is normalized) to those near I~ (relative 
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3.2. Holonomy 

In this subsection, we derive one of the two main equations for the red shift. 
This equation arises from a simple holonomy argument, and at first seems to be 
all that is necessary for understanding the red-shift. We shall see, however, in the 
next subsection, that there is a subtlety, and this equation must be supplemented by 
another in order to understand the black-hole state. 

The red-shift arises by comparing the result of propagating light rays through 
space-time from to with the parallelism at infinity. Mathematically, this 
corresponds to a holonomy, that is, parallel propagation around a certain closed path 
Tu- In fact, it will be convenient to begin by defining the oppositely oriented path 
— r„, as follows. 




the path at infinity 



Figure 2. The path r„ for 
defining the holonomy is the geodesic 7„ 'closed by a path at infinity'. 



to which V is normalized.) 
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(a) Fix some angular coordinates (6'o,</?o): and for each u foUow the radial 
null geodesic 7u(s) inwards from (u, Oq, (po) through the space-time and out 
to {v,9i,(fii) on I~ , where v = v{u), with v{u) the mapping of surfaces of 
constant phase. 

(b) Close the path "at infinity." What is really meant by this is that we use 
the parallelism at infinity to identify the tangent space at {v,9i,<f>i) with that 
at (w, 6*0, v^o)- We abuse terminology slightly by imagining that this parallelism 
arises from transport along some mathematically fictitious path at infinity. 

Transport and parallel along 7u(s). 

Let H°'i,{u) be the holonomy around the path r„, or rather the part of the holonomy 
acting on vectors orthogonal to the surfaces of spherical symmetry, that is, those 
spanned by 1°- and n". Since parallel transport must preserve null vectors, the vectors 
and n"" must return to multiples of themselves, or of each other. The latter possibility 
would require a non-trivial topology and will be ignored. Thus we must have 

= /i^nb + h-^n% (2) 

for a scalar function h{u) which characterizes the holonomy. 

We now derive a differential equation for h. The relative increment in the holonomy 
as we move from r„ to Vu+su is 5u times 

H-^H = HH-^ = (h/h) {rrib - n%) . (3) 

On the other hand, this quantity may be computed by integrating the curvature tensor 
along 7„. If we parallel transport Z" and along 7„, then l"- will be tangent to 7^ and 
n^-Su will be a connecting vector to ju+Su, so we shall have 

H-^H = J Rpgh'^Fn'^ ds (4) 

where s is an affine parameter normalized to PVaS = 1 (and we understand that the 
integration implicitly relies on parallel transport to identify tensor fields along 7^). 
Comparing (3) and (4), we have 

h/h= I i?p^6«Z%%„ZMs . (5) 

We next recast this as a differential equation for v{u). 

By definition, the result of transporting the vector = du at X"*" backwards along 
7„ to X~ is a vector v{u)do at X~ . If this transport is then closed "at infinity," we 
come back to v{u)du, and thus 

t;(w)n" = H-^n" = hn'' , (6) 

or 



v{u) = h . 



(7) 
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Note that this means v{u) > for all u. Thus we have 

(log [ Rods, (8) 

where Rd = Rabcdl°"n^l'^n'^ is the sectional curvature. This is the first difi^erential 
equation we have been seeking. 

3.3. Interpretation and Geometry: Evolution of ds 

The equation (8) evidently governs the evolution of the red-shift, and so at first 
one would think that this is all that is necessary for understanding the approach to 
the black-hole state. However, as mentioned above, there is a subtlety. To understand 
this, note the following: 

(a) The sectional curvature Rd is independent of the choices of scale of 
and as long as the normalization laU"' = 1 is preserved. Thus Rd is well- 
defined by the local geometry of the space-time, and independent of issues of 
normalization at infinity. 

(b) On the other hand, the measure ds is normalized by the requirement 
(""VaS = 1. In other words, as one-forms on the geodesic 7^, we have 

ds = Ha at X+. Notice this means that ds is not normalized relative to the 
Bondi frame at X~, rather it has been red-shifted. 

(c) The measure ds is dynamic; it evolves as u increases. 

The key feature of the geometry turns out to be (c). We shall see that as the 
geodesic 7^ approaches the event horizon, the measure ds tends exponentially to zero 
except on the portion of 7„ connecting T+ to the horizon. It is this, coupled with 
equation (8), which gives rise to the universality of the collapse, for it suppresses all 
contributions of the geometry from times prior to the formation of the hole. 

We may work out the dynamics of ds easily enough. Let us put 

n^V^r = ci^ 

(9) 

(Here q, c„ would be written 23^7, in the Newman-Penrose formalism. However, 
we shall reserve 7 for the geodesic.) Then short computations give 

rv„Q = = -rv„c„ , (10) 

and so 

Cliluis)) = -CnMs)) = - / Rd{1u{s')) ds' . (11) 

The evolution of ds is given by 



Cnds = n^Vfods = c„ds . 



(12) 
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Equation (12) shows that as long as c„ is negative and bounded away from zero 
in a neighborhood of the event horizon, the affine measure ds along the null geodesies 
will tend exponentially quickly to zero as the event horizon is approached. Likewise, 
the measure ds on the portions of 7^ tending to the past of the event horizon will 
vanish exponentially quickly, because they are related to that on the event horizon by 
finite parallel-transport factors.^ 

Note that the argument above shows that, as the horizon is approached, we have 

^ I-yu 



4. Analysis and Applications 

In this section, we look at the consequences of the formulas derived above, and in 
particular study what happens when a black hole forms. We first summarize the main 
definitions and formulas. 

The function v{u) is the mapping of surfaces of constant phase, so v{u) is the 
red-shift of a wave packet propagating through the space-time along the radial null 
geodesic 7^ terminating at X+ at retarded time u. The evolution of the red-shift is 
governed by equation (8) 

{logv{u)y=v/v= / Rods. (8) 

Here Ru = Rabcdf-^'n^l'^'n''' is the sectional curvature, the null vectors (tangent to the 
geodesic) and normalized by /an" = 1. The definition of Rd does not depend on 
further specializations of 1°" and n"", but other quantities in general do, and we assume 
that 1°' and n" are transported parallel along the geodesic and are normalized at X"*". 
The measure ds is the restriction of the one-form Ua to the geodesic, and so it is also 
transported parallel along the geodesic and normalized at T"*". The evolution of ds is 
given by equation (12) 

Cnds = n^Vbds = c„ds , (12) 



where 



Cn{lu{s))= RD{lu{s'))ds' . (11(b)) 



There is implicit here an assumption that the geometry becomes asympotically flat at a reasonable 
rate in a neighborhood of I~ , so that the exponential decrease of ds is not compensated by a growth 
of curvature in a neighborhood of T~ . One would certainly expect this assumption to hold for a 
space-time corresponding to an isolated collapsing object. 
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4-1- Incoming and Outgoing Contributions 

Our equation (8) for the evolution of the red-shift involves an integral over the 
whole of the geodesic 7^. Near an event horizon, the contributions from different 
portions of this are very unequal. Essentially, as u increases, it is the portions of the 
geodesic connecting the outside of the region where the black hole will form to 
which are important, the others being suppressed in the limit. (This portion may 
be thought of as corresponding to Stage I of the discussion in Section 2, whereas the 
contribution of the earlier part of the geodesic corresponds to Stage II.) 

It is therefore convenient to re-write equation (8) by dividing the integral into 
two parts. We choose a point 7u(s(t(.)) on the geodesic 7^, and divide the geodesic 
into a ray 7u_|_ outwards from 7u(s(ti)) to X+ and a ray 7^- from I~ to 7„(s(u)). (We 
shall discuss the choice of 7u(s(w)) shortly.) It will also be convenient to normalize 
the measure on 7„_ relative to I~; this measure is ds_ = {v)'~^ds. We write ds+ = ds 
for the measure on Then we let 

a±(u) = / RDds±. (13) 

It should be emphasized that as long as there is no extreme injection of matter from 
T~ , or ejection of matter to X"*", one would expect the quantities a±{u) to remain finite 
and bounded and (as long as 7u(s(w)) does not go off to infinity) generically non-zero. 
Then equation (8) becomes 

{logv{u))' = a_{u)v + a^{u) (14) 

or equivalently 

V — a_{u){vy + a^{u)v . (15) 



4.2. Phase Portrait 

Equation (15) can be regarded as an evolution equation for v. If we temporarily 
ignore the w-dependence of a±, then the equation has two fixed points, one at = 
and the other at v = —a^/a^. The fixed point v = is attractive for a+ < and 
repulsive for a+ > 0; the other fixed point always has the opposite character. 

In general, there are no restrictions on the signs of a±. The stationary case is 
exceptional in that it turns out that there a_ = a+ = 0, if 7„(s(tt)) is chosen to be the 
spatial origin.^ This is clearly a non-representative situation as far as the behavior 
governed by equation (15); this points up the importance of studying the generic, 
non-stationary case. 

^ Here is an outline of the proof. If there is a timeUke Kilhng vector then the world-hne of the 

spatial origin must be a geodesic with tangent so ^''^bC" vanishes at the origin. However, a short 
computation by resolving into components along 1°" and n° and using Killing's equation shows that 
^*Vj^" comes out to be c„ = times a non-zero spacelike vector. Thus = 0; by symmetry, then 
a_ = 0. 
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According to the sign of a+ and the position of v relative to —a^/a_, then, various 
sorts of quahtative behavior are possible. (Remember that v > always.) 
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dv/du Fig- 
ure 3. Some possible instantaneous phase portraits for v. The fixed points at zero 
and —aj^/a_ are indicated, with the arrows near them indicating their attractive or 
repulsive character. Cases (a) and (d) correspond to collapse, with v being driven to 
zero. Case (b) is stable, with v tending to track —a_^/a_. In case (c), one would have 
V driven towards +oo. 



Approach to an event horizon means v — > 0, and in the simplest cases one would 
expect this to mean that a+(tt) becomes negative and that va+ dominates (?))^a_. 
Precisely, if there are constants o and k such that a+(w) < a < and (if a-{u) > 0) 
\va^/a+\ < k < 1, then we shall have v/v < a(l — k), and hence i) ^ at least as fast 
as ~ expa(l — k)u. One can think of v as collapsing to zero in this case. (See figure 
3(a),(d).) 

On the other hand, if a+ > and a_ < 0, then v = is a repulsive fixed point and 
V = —a+/a- is an attractive one. This is a stable situation. No collapse will occur, 
and V will tend to track the value —a+/a^. (Figure 3(b).) 

A situation with a+,a_ > would drive v larger; if this were to persist it could 
drive v +oo in finite retarded time. Indeed, by comparing equation (15) with 

v^{vfa^{u), (16) 

which is easily integrated, one sees that v +oo at or before the time u such that 

a_{u')du' = l/v{uo) . (17) 



pu 

J Uo 



It is hard to see that such a situation could be compatible with physical requirements 
like global hyperbolicity and reasonable asymptotic behavior at X. It would mean 
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that by the retarded time u, distant observers have aheady received all of the data 
incoming from I~. (Figure 3(c).) 

One can also have more complicated situations, in which the signs of a± change. 
Some qualitative understanding of these can be pieced together from the observations 
just given, but in general one needs to understand the functions a±{u) quantitatively. 
For example, we saw in the previous paragraph that i) could be driven to +00 in finite 
time. This would mean that whatever a±{u) would do beyond that would be "too 
late" to affect the behavior of it. 



4-3. Approach to An Event Horizon 

We have already indicated that in the simplest cases one would expect approach 
to an event horizon to correspond to a situation where a+(ti) < a < (and 
\va^/a+\ < /c < 1 if a_ > 0). Here we discuss this situation in more detail. 

(Wc first mention that one can verify that these inequalities hold in a simple 
example. This is constructed as follows. One starts with a Schwarzschild solution 
of mass m_ and static, say perfect fiuid, interior. Then one introduces an incoming 
pulse of radiation of energy Am at a fixed advanced time (say f = 0), so that one 
has in effect a time-reversed Vaidya solution in the exterior of the fluid. If one takes 
7„(s('u)) to be the spatial origin, then one finds 

a^{u) = 

(18) 

a4n) = -^ + ^;— 

where m+ = m_ + Am is the total mass and ro is the value of the radius at which 
7_ crosses v = 0. We have a+ < 0. If a black hole forms, then ro — > 2m+ and 
a+ — > — l/4m+, (minus) the surface gravity.) 

Notice that as long as 7„(s('u)) is chosen to continuously approach a finite point 
in space-time as tt — > +00, the integral a_(tt) = / Rnds- will remain bounded.^ 
Since this integral is multiplied by ii in equation (14), it will not contribute to the 
limiting behavior of v/i). This is true not just for 7„(s(w)) being, for example, the 
spatial origin, but even if 7„(s(w)) tends to any finite point on the event horizon. This 
shows that the limiting behavior of v/i) is independent of all details of formation of 
the black hole, and even of the geometry of any finite point on the horizon. 

While this behavior is easy enough to understand analytically, it is not captured 
very well by any of the standard infinite regimes common in relativity. What is relevant 
here is compressed to a point in the usual conformal diagram: it is the "gap" between 
the event horizon and I'^. 



^ Again, this is true with mild asymptotic assumptions at X~ . For example, assuming that Rj^ has the 
usual peeling behavior uniformly in v for a neighborhood of I~ near the limiting value of v would be 
sufficient. 
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Figure 4. The regime rele- 
vant to the asymptotic behavior of the fractional acceleration of the red-shift v/v-, 
indicated by arrows. It is the gap between the event horizon and future null infinity. 
This regime is compressed to a single point in a conformal diagram, such as this one. 



Let us recall that the evolution of ds is given by 



£„ds = ■nJ'Vbds = c„ds 



where 



c„(7„(s)) 



1+ 



RDhu{s'))ds'. 



(12) 



:il(b)) 



In the situation we are considering, this approaches a constant negative value if 
7u(s) approaches a finite point on the event horizon. This means that ds vanishes 
exponentially quickly as the event horizon is approached along the vector field n°'J 



^ Since n° is itself being exponentially compressed by the same factor, the measure ds would vanish 
linearly in local coordinates. 
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This is really the feature which makes the limiting behavior of 

v/v= / Rods (19) 



independent of all details of the formation of the hole or of the geometry of finite 
portions of the horizon. 

Let us recall that this subsection has been concerned with the formation of an 
event horizon only in the simple case a+(u) < a < 0, and this probably corresponds to 
what happens in a "normal" collapse of a spherically symmetric star. However, it is 
quite possible that other, more complicated, stituations occur, for example in critical 
collapse. 

4 ■4- Surface Gravity 

Usually, the surface gravity of a black hole is only considered to be defined if 
the hole is stationary. Indeed, the surface gravity is defined in terms of the Killing 
vector associated with the stationarity. We shall show here that in the case of a 
stationary black hole, the surface gravity is simply (minus) v/v, and so this quantity 
is a candidate for the notion of "surface gravity" in general. 

Let us write 

^« = Ar + un"" , (20) 

for scalar functions A and u. (These are not the Newman-Penrose quantities often 
denoted by those symbols.) We have A, z/ > in the stationary region, and A, — > 1 
at X"*" (normalization). The components of Killing's equation in the n" plane are 

rvajv = 

qA + n"VaA = (21) 

rV„A + Cniy + n^Vai^ = 

The first implies u = 1 everywhere. (The remaining system is of course 
over deter mined, giving a consistency requirement for a Killing vector to exist.) Now, 
we may use these equations in any of the standard formulas for the surface gravity. 
For example, one has 

- (l/2)V„e'6 = (22) 

as the horizon is approached. Contracting this with n", say, and substituting (21) 
gives 

K — ci = —Cn at the horizon. (23) 
We saw in the previous subsection that this limiting behavior of c„ was 

v/v= Rods, (19) 
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and thus we identify (minus) this with the surface gravity, whether the exterior is 
asymptotically stationary or not. 

Two remarks are in order. First, in the non-stationary case, one would not usually 
expect —v/i) to tend to a constant value, and this is the reason for the terminology 
"limiting behavior" rather than "limiting value." Second, it is arguable whether this 
quantity ought to be called the surface gravity, since it depends on the global structure 
and not just on the local geometry of the horizon. (However, this is true even when the 
surface gravity is defined by means of a Killing field, for then the definition depends 
on the normalization of the field at infinity.) In any event, the identification made 
here of this quantity with —v/v shows that it is interpretable as a relative acceleration 
even in the non-stationary case; whatever one calls it, is it the key quantity in the 
geometric-optics approximation to scattering by the black hole. 

5. Comments on the Non— Spherically Symmetric Case 

Perhaps surprisingly, much of the analysis above can be extended to the non- 
spherically symmetric case without great difficulty. The results arc formulas for 
the gravitational acceleration of wave-vectors by the black hole; these are no longer 
quantifiable by a scalar red-shift, however. We shall only sketch these ideas here. 

Suppose that (M, gab) is a general space-time, which is well-enough behaved that 
I~ and X+ exist as null hypersurfaces, and 7 is a null geodesic from X~ to 
(Actually, all considerations here are local to the geodesic, so it is enough to have 
small portions of X"^ near the endpoints of the geodesic; one does not need much of 
the global structure of I.) 

Let be a parallel propagated tangent vector to 7„. This will be interpreted as 
the nominal wave vector of a wave packet in the geometric-optics approximation, so 
we want to understand how /" is affected by its passage from X~ to X"^. In order to do 
this, we must have some way of comparing vectors on X~ with those on X"*". In general, 
there is no simple way to do this, for in the presence of gravitational radiation one 
does not have covariantly constant vector fields at X, and this leaves aside the issue 
of relating vectors at X~ to those at X+. 

It is possible to overcome these difficulties by various technical means, but at the 
present level of discussion it turns out that we do not have to address this issue. We 
suppose that a Lorentz transformation L„ is chosen identifying the tangent space at 
the past end-point of 7„ with that at the future end-point. In general, there will be 
different choices of L„ possible, and presumably the correct choice in any situation is 
dictated by the physics. 

We may then define a closed path as before, one portion of which is the geodesic 
7„, and the other portion of which is a mathematical fiction, parallel transport along 
the second portion being effected by L~^. Let us write P{u) for the parallel transport 
along 7^; then the holonomy around the closed path beginning and ending at the 
end-point at X~ is Q = L'^P. A wave vector If^ at X~ is taken to Z"^^ = P°'hl\a 
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T+. In order to get a measure of the red shift, this must be compared with I/^feifn, the 
result of transporting the vector from I~ to by the fictional "path at infinity." In 
the sense of this fictional transport, and referring all vectors to X+, we may identify a 
"red-shift operator" Q = PL~^. Then we have 

Q = PL-^ - PL-^LL^^ 
= PP-^Q - QLL-^ 

In this formula, we have 

P\{p-'yb= [ Rp,b''l''n''ds (25) 

the integration of the tensor field being done by identifying the tensor field with a 
tensor in the tangent space to the end-point of 7„ on X+ by parallel transport along 

If p is a point at which the event horizon is smooth, and near p one varies the 
geodesic 7„ to approach (from the past) the generator of the event horizon through p, 
then the normal field to the geodesic becomes compressed as p is approached. This 
means that n^, which restricts to ds on the geodesic, tends to zero. Thus contributions 
to PP"^ from the geometry in the neighborhood of the event horizon are suppressed, 
and the late-time behavior of PP~^ depends only on the portion of the space-time 
connecting the event horizon to X"*", in the same sense as for radial symmetry. 

6. The Hawking Process 

The foregoing results have implications for the usual view of the Hawking process. 
These will be outlined here, although there are a number of technicalities and caveats 
which we shall only mention, their full treatment requiring more extensive discussion 
than can be given in this space. We consider only the spherically symmetric case here. 

The first point to be made is that to leading order, Hawking's model of the emission 
of radiation by an object collapsing to form a black hole is isomorphic to a moving- 
mirror model in two-dimensional Minkowski space. This comes about as follows. 

The physical question at hand is: given a massless quantum field quiescent in the 
far past and propagating through the space-time of a gravitationally collapsing body, 
what does the field look like at late retarded times far away from the body? This is a 
scattering problem, the object being to work out the field operators near X"*" in terms 
of those near I~. 

We shall restrict our attention to the s-wave sector, which can be shown to be 
by far the main contribution. This leaves us with s-waves propagating through the 
collapsing object. For these, for those field modes of moderate wavelength near T+, the 
geometric-optics approximation is a good one, since these arise from very blue-shifted 
modes near the surface of the collapsing object and on the rest of the space-time. If, 



(24) 
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as usual, we rescale the fields to have finite (operator-valued) limits 0o near I by- 
putting (f) = (po/r, then we have 

4>t{u)^-<l>l{v{u)), (26) 

where v{u) is the mapping of surfaces of constant phase, and the minus sign comes 
from reflection through the spatial origin. 

The relation (26) is formally identical to that for scattering of a massless field on 
one side of a mirror with trajectory v = v{u) (in standard null coordinates u — t — x, 
V = t + x) in two-dimensional Minkowski space. To check that the quantum field 
theories are in fact the same, one has to verify that canonical commutation relations 
and field representations are the same; this is routine, using the fact that the original 
space-time was asymptotically flat near I~ . 

We can immediately write down the expectation of the stress-energy near X"*" from 
knowledge of the corresponding result in the moving-mirror theory (cf. Birrell and 
Davies 1982): 



3 (vV _lv^ 
4:\v) 2 V 



Ub (27) 



for the s-wave contribution to the stress-energy near X+. Of course, this equation 
is only valid insofar as: (a) we may neglect all contributions other than the s-wave 
one; and (b) the geometric-optics approximation (here expressed by equation (26)) 
is legitimate. The usual view is that these do give a good treatment of the leading 
physics. See Page (1976a, 1976b, 1977) for the validity of the s-wave approximation. 
See Hawking (1975, pp. 210-211) for an argument that the inclusion of backscattering 
only alters the result (27) by a geometric factor independent of the details of the 
formation of the hole; for another approach to this issue, see Fredenhagen and Haag 
(1990). 

If there is an event horizon with surface gravity k = —v/i), then at late retarded 
times (and for r ^ m, the mass of the object), we will have 



(T^) = (127rr2)- 



1 2 1- 

-K H K 



Ub. (28) 



In particular, if k — > 0, or even if the time over which k changes signiflcantly is much 
smaller that k~^, we have 

(fr)^(487rr^)-'«'U, (29) 

which is the standard result for asymptotically stationary black holes. But the formula 
(28) is valid even for non-stationary black holes, insofar as the general framework of 
Hawking 's analysis is valid. 

Several comments about this are in order: 

(a) Equation (28) can be viewed as giving us corrections to the adiabatic 
approximation for the variation of the emission rate with the surface gravity. One 
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would expect that approximation to break down when the time scale over which k, 
changes to become comparable to n"^, and that is just what appears in (28). 

(b) In particular, there would be no significant modification of Hawking's 
prediction in the case of an isolated Schwarzschild black hole, since in this case \k\ -C 
until the hole is of Planck dimensions, when any scmiclassical treatment is dubious. 

(c) It is possible for k to be negative; indeed, one would expect k < for an 
accreting black hole. 

(d) It is possible that in some cases (l/4)«;^ + (l/2)«; < 0. In this case would 
correspond to a flux of null negative-energy particles. Such behavior is familiar, if not 
wholly understood, in the moving-mirror case. 

(e) Whatever the sign of k, equation (28) shows that Hawking's prediction will be 
substantially modified if the surface gravity changes very rapidly. 

(f) We can understand the validity of the Unruh (1976) vacuum for computing 
the Hawking effect, in terms of the universality discussed above. The point is that the 
Hawking effect arises from the singular behavior of v, and changing the early part of 
the space-time in any way which only prepends a non-singular transmission of null 
rays will not affect this. Unruh's choice of vacuum is simply such a choice which is 
mathematically natural. 

(g) Finally, one of the most important lessons of the moving-mirror models is 
that it is not possible to understand their energy budget without taking into account 
the energy of the mirror and its driving engine as quantum operators, nor without 
taking into account the entanglement of the field state with that of the mirror and 
driving engine (Parentani 1996, Hclfcr 2001). The entanglement occurs becaiise the 
field energies depend on the mirror and its driving engine, which are themselves 
quantum systems: the field energies are functionals of the acceleration of the mirror, 
which is a certain quantum operator. This quantum character cannot be ignored for 
understanding the system's energetics, because observation of the field energy requires 
driving the mirror/engine system into a superposition of energy eigenstates. If there 
is a parallel with black holes, then the treatment of the collapsing object and its 
space-time geometry as classical is questionable. 
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